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Theoretical  Studies  of  Kinetic  Mechanisms 
of  Negative  Ion  Formation  in  Plasmas 


ABSTRACT 

This  technical  program  constitutes  a  theoretical  research  investigation  of 
the  kinetic  mechanisms  of  negative  ion  formation  in  plasmas.  This  study  was 
directed  toward  elucidating  the  mechanisms  of  the  most  important  volume  dependent 
reactions  that  occur  in  hydrogen-ion,  H"(D"),  source  devices,  primarily  of  the 
Belchenko-Dimov-Dudnikov  (BDD)  type  and  toward  evaluating  other  light  negative 
anions,  such  as  Li",  as  possible  sources.  A  primary  goal  of  this  research 
program  was  to  identify  the  most  important  reactions  leading  to  negative  ion 
production  or  destruction  and  to  estimate  these  reaction  cross  sections  and 
product  branching  ratios.  A  further  goal  was  to  explore  new  chemical  sources  for 
the  production  of  light  mass  negative  atomic  ions.  The  results  of  this  program 
furnish  data  and  provide  direction  for  more  detailed  investigations  into  the 
kinetics  of  both  gas  phase  and  gas-surface  reaction  rates  of  importance  in  ion 
source  devices  and  provide  input  for  reliable  modeling  of  such  systems.  This 
investigation  was  carried  out  using  quantum  mechanical  methods.  Both  ab  initio 
and  density  functional  approaches  were  employed  in  these  studies. 
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STATEMENT  OF  WORK 


The  contractor  shall  conduct  a  theoretical  research  investigation  of  the 
electronic  structure  of  negative  molecular  ions  and  the  kinetic  mechanisms  of 
negative  ion  formation  in  plasmas.  This  investigation  shall  include  the 
calculation  of  potential  energy  hypersurfaces  for  ion-atom  and  ion-molecule  gas 
phase  interactions  that  occur  in  magnetron-type  negative  ion  source  devices.  The 
contractor  will  perform  theoretical  studies  of  H  surface  adsorption  to  reduce  the 
surface  conversion  mechanism  for  H  production  into  a  series  of  well 
characterized  elementary  processes.  In  addition,  calculations  shall  be  performed 
to  estimate  the  rates  of  the  most  important  volume  dependent  reactions  leading  to 
H~  production  or  destruction  in  such  devices.  Experimental  data  collected  at 
LBL,  LASL  and  AVCO-Everett  from  their  ongoing  studies  with  negative  ion  source 
devices  shall  be  iteratively  examined  and  integrated  with  the  results  generated 
from  this  theoretical  investigation.  The  computations  shall  be  carried  out  with 
ab  initio  and  density  functional  methods  using  electronic  structure  and  kinetic 
cross-sec t ion  codes  that  have  previously  been  developed  or  modified  at  this 
Center.  Specifically,  this  investigation  shall  comprise  the  following: 

a.  Quantum  mechanical  calculations  of  the  kinetic  cross-sections  for  volume 

dependent  reactions  leading  to  the  collisional  formation  and  loss  (destruction)  , 

of  h".  ; 

i 

i 

b.  Calculation  of  dissociative  attachment  cross  sections  for  e  +  Li2  i 

-*•  Li  +  Li  as  a  function  of  collisional  energy  and  rovibrational  energy  of  the 

Li 2  molecule .  j 

c.  Development  of  a  computer  code  and  calculation  of  the  electronic  ] 

excitation  cross  sections  for  Li2  via  electron  collisions. 

d.  Continued  studies  of  negative  ion  formation  on  surfaces  using  large  [ 

cluster  representations  with  fractional  hydrogen  coverage.  J 
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BACKGROUND 


Recent  successes  in  the  Tokamak  program  and  in  other  controlled  thermo¬ 
nuclear  research  programs  (Ref.  1)  have  focused  attention  on  the  problem  of 
developing  an  efficient  high-energy  particle  beam  source.  For  several  applica¬ 
tions,  a  neutral  beam  with  energies  above  200  keV  is  desired.  The  acceleration 
of  negative  ions  (primarily  H-  or  D")  to  such  energies,  followed  by  neutraliza¬ 
tion  through  a  stripping  reaction,  appears  at  the  present  time  to  be  the  most 
efficient  approach  for  producing  a  high  energy  neutral  particle  beam. 

In  another  area,  the  feasibility  of  particle  beams  (both  charged  and 
neutral)  as  military  weapons  has  been  under  study.  The  proposed  endoatmospheric 
applications  require  high  intensity  sources  and  mainly  focus  on  the  problems  of 
beam  stability  and  propagation  characteristics.  Exoatmospheric  applications 
require  lower  intensity  sources  but  of  very  high  quality.  Design  goals  are 
highly  collimated  beams  with  a  narrow  energy  spread. 

A  magnetron-type  negative  H"  source  has  been  reported  by  Belchenko,  Dimov 
and  Dudnikov  (BDD)  (Ref.  2)  that  has  produced  H"  current  densities  of  several  A 
ctn“2.  This  device  operates  as  a  plasma  discharge  in  an  atmosphere  of  hydrogen 
gas  with  cesium  or  other  alkalis  present  at  ~  0.01  percent.  The  mechanism  for 
the  production  of  H"  is  believed  to  involve  a  catalyzed  surface  reaction  whereby 
H"  ions  are  produced  by  transfer  of  electrons  from  Cs  atoms  that  are  adsorbed  on 
the  cathode  surface  (Refs.  3,  4).  The  detailed  kinetic  mechanisms  of  such 
surface  reactions  are  still  uncertain  and  parametric  experimental  studies  are 
currently  underway  at  LASL,  LBL,  IRT  and  at  other  laboratories  to  elucidate  the 
mechanisms  and  operating  characteristics  of  BDD  and  similar  devices  and  to 
develop  information  for  their  scale-up  to  higher  current  densities.  Diagnostics 
of  H"  source  devices  are  also  underway  at  Brookhaven  (Ref.  5)  using  beam  probe 
and  spectroscopic  techniques. 

Concurrent  with  these  surface-plasma  reactions  are  several  electron-H 
volume-dependent  processes  that  may  lead  either  to  the  production  of  H"  or,  in 
reverse,  may  act  as  important  destructive  processes  of  the  negative  ions  after 
they  are  formed.  Photodetachment  experiments  by  Bacal  and  Hamilton  (Ref.  6)  in 
hydrogen  plasmas  indicate  H"  densities  100  times  larger  than  that  predicted  from 
simple  electron  attachment  mechanisms.  Further,  these  experiments  indicate  a 
nonlinear  dependence  of  the  production  of  II"  on  electron  density,  at  least  for 
densities  less  than  lO1^  cm"3.  Several  mechanisms  (Ref.  7)  have  been  proposed  to 
explain  these  volume-dependent  H"  production  processes  but  the  cross-sections  for 
many  of  the  reactions  leading  to  H"  are  uncertain. 
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By  a  large  margin,  the  major  thrust  of  the  development  of  particle  beam 
technology  has  been  experimental  in  nature.  Theoretical  studies  of  surface 
catalyzed  reactions  leading  to  H"  production  have  been  undertaken  by  LLL 
(Ref.  4),  but  only  a  relatively  small  theoretical  effort  (Ref.  8)  has  been 
directed  toward  detailed  calculations  of  cross-sections  for  volume-dependent 
reactions.  Because  of  the  difficulty  of  conducting  experiments  to  measure  these 
unknown  cross-sections,  the  development  of  a  sound  theoretical  method  for 
calculating  these  state-to-state  processes  appeared  desirable.  Although 
relatively  little  work  of  this  nature  has  been  done  in  the  past,  enough 
theoretical  work  is  available  to  indicate  that  the  development  of  such  procedures 
can  be  made  practical,  particularly  if  good  wave  functions  and  potential  energy 
surfaces  are  available  for  the  interacting  species. 

The  original  goal  of  this  technical  program  was  to  delineate  the  most 
important  volume  dependent  reactions  that  lead  to  H-(D~)  production  (or 
destruction)  and  to  calculate  the  energetics  and  cross-sections  of  such  reactions 
as  a  function  of  system  parameters.  Attention  was  also  directed  toward  those 
kinetic  processes  that  are  of  importance  in  high  pressure  ion  source  devices, 
particularly  of  the  BDD  type.  It  was  intended  that  this  theoretical  program 
should  compliment  experimental  studies  with  negative  ion  source  devices.  A 
cooperative  program  of  ongoing  experimental  research  and  our  theoretical  studies 
are  of  mutual  benefit  and  result  in  a  better  understanding  ot  the  fundamental 
physical  processes  occurring  in  ion  source  devices. 

A  further  goal  of  this  research  program  is  to  explore  new  chemical  sources 
for  the  production  of  light  negative  atomic  ions,  in  particular  H  and  Li  . 
Current  studies  include  the  formation  of  H  or  D  via  photodissociation  of 
unsaturated  hydrocarbon  negative  ions  and  the  formation  of  Li  via  dissociative- 
attachment  of  electrons  to  molecular  Li2-  This  latter  process  appears  feasible 
based  on  our  current  knowledge  of  the  energetics  of  this  system.  Our  recent 
research  efforts  have  been  directed  toward  three  general  areas:  mutual 
neutralization  and  collisional  dissociation  reactions  in  negative-ion  plasmas, 
analysis  of  dissociative-attachment  in  the  e  +  Li 2  system,  including  the 
calculation  of  electron  excitation  cross  sections  for  e  +  Li2,  anci  further 
studies  of  surface  catalyzed  reactions  for  the  production  of  H  .  The  negative- 
ion  plasma  studies  entail  analysis  of  the  kinetics  of  volume  reactions  leading  to 
the  destruction  of,  for  example,  H  and  Li  .  The  e  +  Li 2  dissociative  attachment 
studies  parallel  the  work  of  Bardsley  and  Wadehra  (Refs.  8,9)  on  the  e  + 
reaction  and  are  being  carried  out  in  collaboration  with  Professor  Wadehra. 

In  addition  to  studies  of  thermodynamically  stable  negative  i  >ns ,  there  has 
b  ^n  recent  interest  in  the  metastable  or  resonant  states  in  these  systems. 
Compound  states  of  atomic  negative  ions  have  been  a  subject  of  interest  in  atomic 
physics  for  over  forty  years.  The  first  observation  of  a  compound  negative  ion 
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state  was  the  detection  of  long-lived  He-  ions  in  a  mass  spectrometer  by  Hiby 
(Ref.  10)  in  1939.  Holoien  and  Midtdal  (Ref.  11)  theoretically  showed  that  He- 
(ls2s2p)4P°  was  stable  relative  to  He(ls2s)^S  and  thereby  offered  an  explanation 
of  Hiby's  observation  of  long-lived  He-  ions.  In  1963,  Schulz  (Ref.  12)  reported 
evidence  for  a  scattering  resonance  in  He  at  ~  19.3  eV  which  was  attributed  to 
He~( ls2s 2) 2S .  This  was  the  first  experimental  observation  of  an  electron  scat¬ 
tering  resonance  in  an  atomic  system.  Since  that  time,  compound  states  of  atoms 
have  been  seen  in  many  elastic  and  inelastic  channels  in  electron-atom  scattering 
(Ref.  13),  in  ion-atom  scattering  cross  sections  (Ref.  14),  in  laser  photodetach¬ 
ment  cross  sections  (Ref.  15),  and  in  dissociative  electron  attachment  studies 
(Ref.  16,  17).  Compound  states  of  doubly  charged  negative  ions  have  also  been 
detected  (Ref.  18). 

The  spectroscopy  of  metastable  negative  ions  has  been  the  subject  of  several 
experimental  and  theoretical  investigations  (Refs.  19-23).  Although  much  atten¬ 
tion  has  been  given  to  the  problem  of  atomic  anion  formation,  more  recent  studies 
have  examined  diatomic  anions,  especially  those  o'  the  lighter  elements 
(Refs.  24-27).  In  a  recent  study  (Ref.  25),  the  electronic  structure  of  the 
newly  discovered  (Ref.  24)  metastable  He£  ion  was  examined  using  Cl  wavefunc- 
tions.  This  study  indicated  that  only  the  state  of  He!  was  metastable  rela- 

tive  to  He2  [ a J  Eu ]  +  e.  More  recently,  by  double  charge  exchange  of  He2  ions  in 
lithium  vapor,  Kvale,  et  al  (Ref.  28)  have  examined  the  autodetachraent  spectrum 
for  ^He^  and  He^  ions  and  have  observed  two  autodetachment  peaks.  The  first 
(low-energy)  peak  was  identified  as  autodetachment  from  the  first  vibrational 
level  of  He 2  [  IT  ]  to  He2  [a'*!*]  +  e  and  the  second  as  autodetachment  to  the  He2 
[XE*i  +  e  repulsive  continuum,  similar  to  that  observed  by  Bae,  et  al 
(Ref!  24). 

The  coupling  of  this  experimental  data  to  theoretical  studies  of  the  elec¬ 
tronic  structure  of  the  He2  anion  has  yielded  a  good  understanding  of  the  spec¬ 
troscopy  of  this  system.  Because  of  the  difficulties  in  conducting  experiments 
to  measure  kinetic  processes  in  these  metastable  and  compound  ion  systems,  the 
development  of  sound  theoretical  methods  for  calculating  the  positions  and  widths 
from  analysis  of  their  electronic  structure  appears  desirable.  In  this  direc¬ 
tion,  we  have  undertaken  a  series  of  calculations  for  the  He/He  /He2  and  Li/Li 
/Li^  systems  in  an  effort  to  understand  their  resonance  structure. 

The  calculations  for  the  He’/He^  system  were  carried  out  to  gain  a  better 
understanding  of  the  experiments  of  Bae,  Coggiola  and  Peterson  (Ref.  24),  who 
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were  the  first  to  report  observations  of  the  raetastable  He2  ion,  and  the  more 
recent  autodetachment  experiments  reported  by  Kvale,  et  al  (Ref.  28)  for  this 
system.  Our  results  show  that  only  the  lowest  of  He2  is  metastable  relative 

to  He7  [a^E  +]  +  e  and  that  the  lowest  lying  ,  ^E*  and  ^E*states  are  all 
autoionizmg .  However  since  the  ^5/2g  sPin*“or^lt  state  has  a  long  400  usee) 
lifetime,  it  should  be  possible  to  laser  induce  transitions  from  this  level 
to  the  upper  resonant  negative  ion  states.  Such  an  experiment,  using  a  tunable 
laser,  is  now  under  development  by  Compton  and  co-workers  at  Oak  Ridge  National 
Laboratory.  A  complementary  experiment  would  involve  the  possible  production  of 
thermodynamically  stable  negative  ions  by  stimulated  radiative  attachment  of 
electrons  to  neutral  atoms  or  molecules.  An  alkali  system  would  appear  to  be  a 
good  choice  for  such  studies  since  both  Li-  and  are  stable  and  exhibit  an 

electronic  structure  which  includes  resonances  that  could  be  reached  by  conven¬ 
tional  tunable  lasers. 


In  a  related  area,  there  is  now  considerable  interest  in  understanding  ion- 
pair  production  in  low  energy  collisions  of  Rydberg  atoms.  Ciocca,  et  al  (Ref. 
29)  and  Cheret  and  Barbier  (Ref.  30)  have  reported  efficient  ion-pair  formation 
in  laser  excited  alkali  systems.  In  contrast,  McGeoch  and  Schlier  (Ref.  31)  find 
no  ion-pair  production,  as  a  primary  process,  in  thermal  energy  collisions 
between  lithium  Rydberg  atoms.  Conventional  curve-crossing  mechanisms  are  not 
capable  of  explaining  these  results  and  new  theoretical  developments  are  needed. 
An  experiment  to  study  ion-pair  formation  in  the  reaction  of  He(3S)  atoms  and 
alkali  Rydberg  atoms  has  been  suggested  by  Peterson  (SRI).  This  system  would  be 
an  ideal  candidate  for  a  collaborative  experimental  and  theoretical  study  of  the 
mechanism  of  ion-pair  production. 


The  goals  of  this  research  program  are  to  identify  the  electronic  and 
resonance  structure  of  atomic  and  molecular  positive  and  negative  ions,  to  study 
the  effects  of  high  laser  fields  on  stimulated  autodetachment  and  radiative  elec¬ 
tron  attachment  in  such  systems,  and  to  carry  out  theoretical  studies  of  atomic 
Rydberg  state  interactions  leading  to  ion-pair  production.  This  theoretical 
research  effort  should  complement  the  experimental  studies  of  ion  formation 
processes  that  are  being  carried  out  by  Peterson  and  co-workers  at  (SRI),  by 
McGeoch  (AVCO)  and  by  several  European  research  groups. 
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RESEARCH  PROGRESS 

During  the  course  of  this  contract,  we  have  undertaken  studies  of  several 
light  negative  ion  systems.  As  a  continuation  of  our  work  on  Li-  formation,  we 
have  suggested  a  new  kinetic  route  for  the  formation  of  quiescent  Li~  through 
Rydberg  atom  collisions.  This  new  mechanism  is  now  under  study  by  M.  McGeoch  at 
AVCO  as  an  alternative  to  dissociative  attachment  of  e  +  Li2  as  a  source  for  Li-. 
In  addition  to  our  Li-  studies  we  have  undertaken  an  analysis  of  the  kinetics  of 
He^  autodetachraent .  This  work  is  being  done  in  collaboration  with  J.  Peterson 
and  Y.  Bae  at  SRI  and  has  resulted  in  several  publications  and  joint  technical 
papers.  Finally  we  have  examined  the  negative  hydrogen  ion  cluster,  H^,  to 
determine  whether  this  system  is  thermodynamically  stable  relative  to  H'  +  H^.  A 
summary  of  our  research  is  given  below.  Six  technical  papers  and  nineteen 
presentations  have  resulted  from  the  research  efforts  under  the  present 
contract . 


The  negative  ions  of  light  atoms  are  currently  being  studied  for  their 
possible  application  in  gaseous  discharges,  fusion  plasmas  and  gas  lasers  (Ref. 
32).  One  source  for  the  volume  production  of  atomic  anions  is  the  process  of 
dissociative  electron  attachment  to  molecules  (Ref.  33).  This  process  has  been 
studied  experimentally  for  low  energy  electron-hydrogen  molecule  collisions  by 
Allen  and  Wong  (Ref.  34).  A  parallel  theoretical  study  has  been  reported  by 
Wadehra  and  Bardsley  (Refs.  8  and  9).  More  recently,  McGeoch  and  Schlier 
(Ref.  35)  have  examined  dissociative  attachment  (DA)  in  electron-lithium  molecule 
collisions  and  have  found  large  DA  rates  for  attachment  to  highly  vibrationally 
excited  Li 2  molecules.  The  effect  of  vibrational  excitation  on  the  DA  rates  has 
been  studied  by  Hiskes  (Ref.  36)  for  e  +  H2  collisions  and  by  Wadehra  and  Michels 
(Ref.  37)  for  the  e  +  Li2  system. 

During  the  past  several  years,  we  have  undertaken  a  series  of  calculations 
of  Li2  in  order  to  assess  the  energetics  of  this  anion.  McGeoch  (Ref.  38)  has 
suggested  that  a  Li-  source  should  be  possible  through  a  series  of  reactions 
similar  to  those  found  for  H2 .  The  initial  step  involves  the  formation  of  the 
Li 2  dimer  and  excitation  to  low-lying  excited  states  via  fast  electron  collisions 
or  photon  pumping.  The  latter  process  is  that  currently  being  investigated  by 
McGeoch.  The  reactions  are: 


(fast)  +  Li 2  [X  ll+]  >  Li*  [A  lZ*.  B  lnu]  +  e 
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or 


Li 2  [X  XIg]  +  hv  ■*>  Li*  [A  lZ+,  B  hj  (2) 

The  A  * and  B  states  both  exhibit  fast  radiative  decay  back  to  the  ground 
X  ^  state  but  the  Fr anck-Condon  transitions  favor  a  shift  in  the  population  of 
Li 2  to  higher  vibrational  states: 


B  1nu]  +  Li 


* 

2 


[X  1E*  (v  >  5)]  +  hv 


(3) 


Our  present  understanding  of  the  energetics  of  the  Li2/Li2  system  is  shown 
in  Figs.  1  and  2.  The  dashed  curves  in  Fig.  1  were  estimates  based  on  H2  and  the 
known  molecular  orbital  structure  of  LiT.  Our  most  recent  ab  initio  calcula- 

•  ,  2^+  , »  , 

tions,  which  indicate  that  the  lowest  Eg  state  of  L12,  is  slightly  attractive 
all  the  way  into  the  curve-crossing  with  the  ground  state  of  Li2,  are  shown  in 
Fig.  2.  This  suggests  that  electron  dissociative-attachment  is  possible  for  this 
system  for  vibrationally  excited  Li2  molecules: 

e  (slow)  +  Li 2  [X  XEg]  -►  Li  [2S]  +  Li“[1S]  (4) 


The  final  potential  energy  curves  governing  reaction  (4)  are  shown  in  Fig. 

2.  There  is  still  considerable  uncertainty  in  the  location  of  the  excited  reaso- 
nant  states  of  Li,,  especially  in  the  auto-ionizing  region  above  Li  +  e.  The 
ground  (X  j  state  is  well  characterized  theoretically  (Refs.  39  and  40)  but 
little  is  known  about  the  location  of  the  low-lying  excited  electronic  states 
which  connect  to  Li  [2P]  +  Li  [*S].  In  addition,  by  analogy  with  H2,  low-lying 
electronic  states  of  Li~  with  MO  configuration  (la2  la2  2o  liT2)[2,^Eg]  may,  in 
fact,  actually  lie  lower  than  some  of  the  states  shown  in  Fig.  1.  Such  states 
dissociate  to  very  high  atom-ion  limits  and  their  general  characteristics  are 
uncertain.  The  curves  shown  in  Fig.  2  represent  our  best  estimates  at  the 
present  time  but  more  definitive  calculations  are  required  to  define  the  curve¬ 
crossing  regions  for  dissociative-attachment. 


The  lowest 


state  of  Li2  has 


the  molecular  orbital  configuration 
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2  2  2 

[lo_  1 au  2au  2a„]  for  large  separations.  At  short  internuclear  separations,  this 


©V*'J&  2+  2  1  • 
state  mixes  with  a  second  bound  zEg  state  arising  from  Li[zP]  +  Li~[1S]  with  a 


molecular  orbital  configuration  [fag  la^  2ag  3og].  This  mixing  at  intermediate 

separations  results  in  a  lowering  of  the  interaction  potential  for  the  lowest 
2  +  , 

Eg  state.  Our  preliminary  calculations,  as  shown  in  Fig.  1,  indicate  that  this 
mixing  is  stronger  than  expected  from  MO  considerations  and  lowers  the  ^E*  poten¬ 
tial  curve  into  a  feasible  region  for  dissociative-attachment  to  occur. 


The  resonance  structure  associated  with  electron  scattering  by  the  hydrogen 
molecule  has  been  the  subject  of  many  theoretical  and  experimental  studies  (Refs. 
41-43).  In  contrast,  electron  scattering  by  Li2  has  received  relatively  little 
attention.  Miller,  et  al.  (Ref.  44)  have  measured  total  cross  sections  for  e  + 
Li2  scattering  in  the  0.5-10.0  eV  region  and  find  a  very  large  (~  500  A^)  and 
rising  total  cross  section  for  the  low-energy  region.  These  results  suggest  that 
some  sort  of  shape  resonance  may  be  found  in  this  system  corresponding  to  a  low- 
lying  state  of  Li2>  Such  resonant  states  have  been  observed  in  electron  impact 
excitation  of  Na2  (Ref.  45). 

In  order  to  clarify  the  electronic  structure  of  the  Li2  anion,  a  series  of 
detailed  ab  initio  calculations  of  the  parent  *Eg  state  of  Li2  and  of  several 
possible  symmetries  of  the  anion  were  undertaken.  The  calculational  framework 
that  was  chosen  was  a  configuration-interaction  (Cl)  expansion  built  from  ortho¬ 
normal  Slater-type  orbitals  (STO's).  An  extended  (14a,  6it,  26)  basis  was  opti¬ 
mized  using  several  methods  of  approach.  Complete  orbital  optimization  was 

2  +  __ 

possible  only  for  Li2  and  the  ground  X  Eu  state  of  Li2  since  the  excited  anion 

states  are  either  autodetaching  at  some  critical  internuclear  separation  or  are 

2  +  — 

pure  resonant  states.  The  lowest  Eg  state  of  Li2,  which  is  the  principle 
channel  for  DA  of  e  +  Li2,  is  var iationally  stable  for  large  internuclear  separa¬ 
tions,  and  a  smooth  extrapolation  of  optimized  orbital  exponents  for  this  state 
was  carried  out  into  the  autodetaching  region.  The  purely  resonant  states  of 
Li  and  Li2  were  treated  using  a  modified  nuclear  charge  hamiltonian,  similar  to 
that  described  by  Nestmann,  et  al .  (Ref.  46).  In  this  method  the  resonant  states 
are  variationally  stabilized  by  an  increased  nuclear  charge  and  an  analytic 
extrapolation  is  carried  out  to  determine  the  real  part  of  the  expectation  value 
of  the  unperturbed  hamiltonian. 


Our  calculated  electron  affinity  for  Li,,  derived  from  our  location  of  the 
X  E  state  of  Li9,  is  0.42  eV,  m  good  agreement  with  previous  theoretical 

r  2  , —  .  . 

studies.  The  A  Eg  state  of  Li2  exhibits  a  complex  potential  energy  for 
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inter-nuclear  separations  smaller  than  3.45  A,  where  a  crossing  with  the  X  ^E* 
state  of  Li2  is  found.  This  resonant  state  of  Li2  is  mainly  of  the  Feshbach  type 
and  has  a  long  lifetime  against  autoionization.  We  find  a  small  width  for  decay 
to  the  ground  state  of  Li,,  T(R)  =  0.01  k(R).  Of  special  interest  is  the  low- 
lying  H.  resonant  state  which  lies  approximately  150  meV  above  the  ground  state 
of  Li2.  This  anion  state  has  a  dominant  MO  configuration  of  (log  lau  2o  liru) , 
indicating  a  relatively  large  capture  width  for  low-energy  electron  scattering, 
and  is  probably  responsible  for  the  large  scattering  cross  sections  observed  by 

Miller,  et  al .  (Ref.  44)  for  low  collisional  energies.  The  character  of  this 

2  —  .  u 
state  is  similar  to  the  ITg  resonant  state  of  N2  in  that  the  anion  state  lies 

lower  than  the  corresponding  neutral  parent  for  a  wide  region  of  internuclear 

separations  (3.5  Z  R  Z  9.0  A). 

Further  theoretical  studies  of  these  interaction  potentials  are  required  to 

define  the  overall  reactive  kinetics  leading  to  Li~  production.  A  study  parallel 

to  that  reported  by  Hiskes  (Ref.  36)  for  H2  has  been  carried  out  to  estimate  the 

overall  cross  sections  for  production  of  vibrat ional ly  excited  Li2>  In  this 

study,  cross-sections  for  vibrational  excitation  (VE)  of  the  X^E*  state  of  Li2 

via  formation  of  the  excited  A  ^E+  or  B  electronic  states  were  reported.  For 

u  1  +  u 

vibrational  excitation  through  the  A  Eu  state,  the  cross  sections  are  nearly 
constant  for  forming  3  <  v"  <  9  via  electron  collisional  excitation.  For  photon 
pumping  (670  nm  <  X  <  700  nm)  of  the  A  ^E*  state,  levels  v"  <  9  are  predominantly 
formed.  Vibrational  excitation  via  the  B  ^IIu  state  has  a  lower  probability. 

In  addition,  theoretical  studies  are  in  progress  to  estimate  the  low-energy  elec¬ 
tron  dissociative  attachment  cross-section  for  e  +  Li*  Li  +  Li  .  Our  initial 
results  are  shown  in  Fig.  3  which  indicate  an  attachment  cross  section  >  10-16 
cm2  for  vibrationally  excited  Li2>  The  enhancement  in  the  attachment  cross 
section  with  vibrational  energy  is  similar  to  that  found  for  the  e  +  H2*  attach¬ 
ment  . 

An  alternative  mechanism  for  negative  ion  formation  was  suggested  by  Lee  and 
Mahan  (Ref.  47),  who  proposed  that  ion-pair  formation: 

X*  +  X  *  X+  +  X-  (5) 

should  occur  for  systems  where  the  energetics  for  ion-pair  formation  are  compet¬ 
itive  with  associative  ionization: 
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The  alkali  metals,  including  lithium,  are  therefore  good  candidates  for  such  a 
study  of  ion-pair  formation.  Ciocca,  et  al .  (Ref.  29)  have  recently  reported  on 
the  formation  of  Na~  ions  in  highly  excited  (n=7-40)  Rydberg  atom  collisions. 

In  a  similar  experiement,  Cheret  and  Barbier  (Ref.  30)  have  reported  efficient 
ion-pair  formation  in  collisions  between  ground  state  and  Rydberg  Rb  atoms  formed 
by  laser  excitation.  A  curve  crossing  mechanism  was  suggested  to  explain  these 
results . 

More  recently,  McGeoch  and  Schlier  (Ref.  31)  have  looked  for  ion-pair  form¬ 
ation  in  collisions  between  two  Rydberg  Li  atoms  (n  =  5-9).  The  most  prominent 
processed  observed  was  associative  ionization: 

.  irk  .  irk  .  +  .  _ . 

Ll  ♦  Ll  ♦  L12  +  e  (7) 

rather  than  ion-pair  production,  which  was  not  observable  at  a  rate  ~  300  times 
less  than  associative  ionization.  McGeoch  and  Schlier  followed  the  time  depend¬ 
ence  for  the  formation  of  Li"  and  find  that  it  follows  the  production  route  by 
associative  ionization  of  e  +  Li2>  These  results,  for  the  first  time,  rule  out 
direct  ion-pair  formation  as  an  important  process,  at  least  in  Li  atom 
collisions,  and  suggest  that  the  earlier  work  of  Ciocca,  et  al .  (Ref.  29)  and 
Cheret  and  Barbier  (Ref.  30)  may  be  subject  to  a  different  interpretation. 

In  order  to  determine  the  possible  paths  for  alkali  excited  state  interac¬ 
tions,  we  have  undertaken  a  series  of  ab  initio  calculations  of  the  electronic 
states  of  hi 2  up  to  the  Li[2S(3s)]  +  Li[2S(3s)]  dissociation  limit.  Since  the 
formation  of  the  ion-pair  Li+(1S]  +  I-i"~  1 1 S ]  can  occur  only  for  and  *£+ 
symmetries,  we  have  restricted  our  studies  to  these  two  representations  (see 
Table  1). 

Theoretical  potential  curves  for  the  and  symmetries  of  Li2  were 

obtained  from  valence  configuration  interaction  (VCI)  calculations.  A  64  func¬ 
tion  Slater  orbital  basis,  containing  a  (2slp)  optimized  valence  basis  augmented 
by  a  (3s3p3d2f)  Rydberg  basis,  was  used  in  this  study.  The  basis  was  transformed 
to  symmetry  orbitals  and  Cl  calculations  were  performed  in  the  space  of  all 

symmetry  adapted  configurations  having  the  lCg  orbital  doubly  occupied.  The 
resulting  Cl  expansions  contained  344  and  331  configurations  for  the  *£  and  ^E 
symmetries,  respectively.  As  the  interactions  in  this  system  are  very  long- 
range,  calculations  were  performed  over  a  range  of  internuclear  separations  from 
5  to  50  bohrs. 

The  results  of  our  theoretical  calculations  are  shown  in  Fig.  4  and  Fig.  5 
for  the  1  £+  and  symmetries,  respectively.  The  ionic  curves  which  dissociate 
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to  Li+Ps]  +  Li"[1S(2s2)]  are  clearly  evident  in  both  Fig.  4  and  Fig.  5,  where 
they  exhibit  a  nearly  diabatic  crossing  behavior  with  the  normal  valence  states 
of  Li 2  at  large  internuclear  separations,  but  show  strong  mixing  at  distances 
<10A.  Ion-pair  formation  may  occur  for  long-range  interactions  of  the  type 

Li  [ 2S(2s) ]  +  Li**[ 2S(n=5) ]  *  Li+  [ 1 S ]  +  ^i~ [ 1 S ]  (8) 

but  clearly  there  are  few  effective  curve  crossing  channels  leading  to  Li+  +  Li”. 
The  curve  crossing  mechanism  suggested  by  Ciocca,  et  al.  (Ref.  29),  thus  appears 
not  to  be  an  important  mechanism  for  ion-pair  production  in  the  Li2  system. 

A  second  mechanism  suggested  by  Ciocca  et  al.  (Ref.  29)  for  alkali  negative 
ion  formation  would  involve  excited  state  negative  ion  formation  via 

Li*[2P(2p)]  +  Li**(n>9)  Li+[1S]  +  Li"  [  1P(2s2p)  ]  (9) 

or  from  electron  attachment  to  a  Rydberg  state  of  Li  followed  by  radiative 
stabilization  to  the  ground  Li“[1S(2s2)]  state.  Both  of  these  processes  can  be 
ruled  out  for  the  Li2  system  since  the  autoionization  lifetime  of  the  Li“ 
[1P(2s2p)]  state  is  known  to  be  very  short  (Ref.  48)  relative  to  radiative 
stabilization  and  higher  autoionizing  states  should  exhibit  even  shorter 
lifetimes . 

As  mentioned  above,  McGeoch  and  Schlier  (Ref.  31)  have  recently  found  that 
the  most  prominent  process  occurring  in  lithium  Rydberg  atom  collisions  is 
associative  ionization: 

Lt+Li+Li^+e  (10) 

This  reaction  has  excess  exothermicity  beyond  the  dissociation  energy  of  Li2+. 
Some  mechanism  for  disposing  this  excess  energy  must  be  developed,  either  by 
radiation  or  by  internal  modes.  It  is  very  unlikely  that  the  excess  energy  is 
released  as  kinetic  energy  since  very  hot  electrons  would  be  formed  which  would 
result  in  secondary  ionization  processes.  An  alternate  possibility  is  that  an 
excited  metastable  state  of  Li2+  is  formed  with  sufficient  lifetime  to  be 
detected  and  mass  analyzed.  Detailed  studies  of  such  mechanisms  are  proposed  as 
part  of  a  new  research  program. 


In  addition  to  our  studies  of  the  L^/L^  system,  calculations  were  carried 
out  on  the  He^  raetastable  anion.  The  existence  of  metastable  negative  ions  has 
been  the  subject  of  several  recent  experimental  and  theoretical  investigations 
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(Refs.  49-51).  Although  much  attention  has  been  given  to  the  problem  of  atomic 
anion  formation,  more  recent  studies  have  examined  diatomic  anions,  especially 
those  of  the  lighter  elements  (Refs.  52-55).  In  a  recent  study  (Ref.  53),  we 
have  examined  the  electronic  structure  of  the  newly  discovered  (Ref.  52) 
metastable  He2  ion  using  Cl  wavefunctions .  This  study  indicated  that  only  the 
4Hg  state  of  He2  was  metastable  relative  to  He2  [a  ^Eu+]  +  e,  with  an  indicated 
electron  affinity  of  0.233  eV  and  a  vibrational  frequency  (w  )  of  1783.0  cm-1.  A 

A  +  __  ,  ®  r 

low-lying  Eg  resonance  state  of  He2  was  also  found  but  it  was  not  possible, 
within  the  error  of  the  calculations,  to  ascertain  whether  or  not  this  state 
might  also  be  metastable. 

More  recently,  by  double  charge  exchange  of  He^  ions  in  lithium  vapor, 

Kvale,  et  al .  (Ref.  56)  have  examined  the  autodetachment  spectrum  for  and 

"^He^  ions  and  have  observed  two  autodetachment  peaks.  The  first  (low-energy) 
peak  was  identified  as  autodetachment  from  the  first  vibrational  level  of  He^ 

[4n  ]  to  He2  [a^E*]  +  e  and  the  second  as  autodetachment  to  the  He2  [ *Eg]  +  e 
repulsive  continuum,  similar  to  that  observed  by  Bae ,  et  al .  (Ref.  52).  This 
experiment  by  Kvale,  et  al .  suggested  that  either  r he  low-energy  autodetachment 
was  from  the  first  excited  vibrational  level  of  the  4H  state  as  postulated,  or 
from  a  second  metastable  electronic  state  (possibly  4E*7  lying  in  this  same  ener¬ 
gy  range.  In  either  case,  our  earlier  results,  which  indicated  that  the  first 
excited  vibrational  level  of  the  4IIg  state  was  non-autodetaching  and  that  the 
lowest  4 state  was  purely  resonant  in  character,  were  in  disagreement  with 
these  experimental  observations. 

In  order  to  clarify  the  spectroscopy  of  He2,  a  re-examination  of  both  the 
a  ^  E^  state  of  He2  and  the  low-lying  symmetries  of  He2  was  carried  out.  Our 
previous  studies  of  this  system  were  based  on  Cl  wavefunctions  constructed  from  a 
(14a,  6ir,  26)  STO-MO  basis.  The  atomic  orbital  parameters  were  taken  from  opti¬ 
mized  calculations  for  the  4P°  state  of  He”.  In  these  new  calculations,  a  com¬ 
plete  non-linear  orbital  optimization  was  carried  out  at  the  equilibrium  internu- 
clcar  separations.  This  non-linear  optimization  was  performed  separately  for 
both  the  neutral  and  negative  ion  states  in  an  attempt  to  quantitatively  charac¬ 
terize  the  small  energy  shifts  arising  from  polarization  effects  in  the  formation 
of  the  anion.  This  approach  permitted  examination  of  smaller  changes  in  the 
charge  distribution,  upon  anion  formation,  than  would  have  been  possible  by  fixed 
basis  set  expansions  with  exponents  chosen  to  insure  freedom  from  linear  depen¬ 
dency  problems. 
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For  studies  of  the  lowest  ^E*  state,  a  modified  nuclear  charge  hamiltonian 
was  employed  similar  to  the  method  described  by  Nestmann  and  Peyerimhoff  (Ref. 
57).  Stabilization  of  this  state  (relative  to  He2  [a  ^E*]  +  e)  was  realized  by 
increasing  the  nuclear  charge  by  ~  0.1%.  A  careful  study  of  the  non-linear 
parameters  as  a  function  of  increased  nuclear  charge  clearly  indicated  that  the 
lowest  ^Eg  state  is  purely  resonant  in  character  and  that  this  state  is 
autoionizmg  relative  to  He2  [a3E*]  +  e.  We  find  that  this  lowest  4Eg  state  lies 
0.074  eV  above  He2  [a3E*]  and  has  a  dominant  MO  structure  of  (lso|lsau  2sou)  . 

A  preliminary  analysis  of  the  complex  part  of  the  potential  indicates  that  this 
state  has  a  large  width  for  decay  to  the  He2  [a  ^E*]  +  e  conti-nuum  and  is  unlike¬ 
ly  to  be  observable  on  the  time  scale  of  the  detection  apparatus  in  the  double 
charge  exchange  experiments.  Further  calculations  on  the  ^E*  and  other  low  lying 
resonances  in  He2  are  required  to  define  the  real  and  complex  behavior  of  these 
states . 

The  lowest  ^11  state  of  He2  was  found  to  be  variat ional ly  stable  for  all 
internuclear  separations  that  were  examined.  Our  final  potential  curve  for  this 
state  is  displayed  in  Fig.  6.  We  find  an  electron  affinity  of  0.182  eV  for  the 
state  of  He2,  a  value  ~  50  meV  smaller  than  in  our  previous  calculations 
which  employed  an  STO  basis  that  was  optimized  for  the  4P°  (Is2s2p)  state  of  He~. 
We  attribute  this  smaller  calculated  electron  affinity  to  our  observation  that 
the  use  of  optimized  basis  functions  resulted  in  a  larger  improvement  in  the 
calculated  energy  of  He2  [a  ^E*]  than  to  the  corresponding  He2[ZfII  ]  anion.  Our 
predicted  electron  affinity,  coupled  with  a  calculated  vibrational  spacing  of 
1702  cm-1,  now  indicates  that  the  first  excited  vibrational  level  of  the 
state  of  He2  is  autoionizing  and  is  probably  responsible  for  the  observed  low 
energy  peak  in  the  electron  detachment  studies  of  Kvale,  et  al .  (Ref.  56).  Their 
experimental  estimate  (0.175  eV)  of  the  electron  affinity  of  the  ^Jlg  state  of 
He2,  relative  to  He2  [a  ^E*],  and  our  calculated  value  (0.182  eV)  are  within  the 
uncertainties  of  both  the  experimental  measurements  and  calculations. 

In  a  related  study  at  SRI,  Bae,  Coggiola  and  Peterson  (Ref.  52)  have 
measured  the  energy  spectrum  resulting  from  electron  autodetachment  from  the 
lowest  4IIg  state  of  He2“.  This  continuum  spectrum  exhibits  a  single  peak  at 
15.7  eV  resulting  from  a  vertical  transition  from  He2“  to  the  repulsive  ground 
state  of  He2<  Both  the  shape  and  position  of  this  autodetachment  spectrum  were 
verified  in  later  studies  reported  by  Kvale,  et  al .  (Ref.  56).  However,  neither 
the  shape  or  position  of  the  experimental  spectrum  seemed  to  fit  predictions 
based  on  our  He2~  calculations  and  the  generally  accepted  ground  state  He2 
potential  reported  by  Foreman,  et  al.  (Ref.  58). 
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In  order  to  determine  the  source  of  this  discrepancy,  we  have  analyzed  the 
intensity  distribution  of  the  Hopfield  helium  continuum  reported  by  Huffman,  et 
al .  (Ref.  60).  This  continuum  emission  arises  from  the  transitions, 

He*  [A 1  ]  -*■  [X1!*]  and  He2  [D1^]  -*■  [X^I*],  where  several  vibrational  levels  of 

the  radiating  states  contribute  owing  to  the  nature  of  the  discharge  formation  of 
the  He2  excimer.  The  spectroscopic  data  for  both  the  A1!*  and  states  of  He2 

are  accurately  known  from  the  work  of  Tanaka  and  Yoshino  (Ref.  61)  and  Ginter 
(Ref.  62),  leaving  any  uncertainty  in  the  analysis  of  the  Hopfield  continuum 
solely  in  the  repulsive  character  of  the  ground  state  potential  of  He2 .  Since 
the  equilibrium  separation  of  the  A*E*  state  of  He2  (1.0406  A)  is  close  to  the 
predicted  value  for  the  '♦H  state  of  He^  (1.064  A),  an  analysis  of  the  He2 
continuum  intensity  should  be  a  sensitive  test  of  the  ground  state  potential  in  a 
region  where  discrepancy  between  theory  and  experiment  has  existed  (Refs.  63-65). 

The  problem  of  calculating  the  bound-free  emission  spectra  for  a  diatomic 
molecule  has  been  previously  treated  in  detail  (Refs.  66  and  67).  For  an 
analysis  of  the  spectral  intensity  distribution  arising  from  the  A*E*  ->•  X*E*  and 

-*■  X1!^  transitions  in  He2j  accurate  numerical  wavefunctions  are  required  for 
both  the  bound  vibronic  levels  of  the  upper  radiating  states  and  the  continuum 
dissociating  state.  Such  wavefunctions  describing  the  nuclear  motion  of  the 
system  can  be  accurately  generated,  provided  reliable  potential  energy  functions 
are  available  for  the  electronic  states  involved  in  the  transition. 

For  the  bound  A1!*  and  states,  accurate  RKR  potential  energy  functions 

were  generated  based  on  the  latest  analysis  of  the  spectroscopic  constants  given 
by  Huber  and  Herzberg  (Ref.  68).  For  the  ground  state  potential  of  He2,  we  have 
used  three  different  functional  forms.  The  first  was  the  exponential  fit  derived 
by  Foreman,  et  al .  (Ref.  69),  based  on  the  their  careful  analysis  of  total 
scattering  cross  sections  for  He-He  collisions.  The  second  potential  employed 
was  the  ab  initio  result  recently  reported  by  Ceperley  and  Partridge  (Ref.  70), 
based  on  quantum  Monte  Carlo  calculations  at  short  internuclear  separations 
(0.5  A  <  R  <,  1.8  A).  This  theoretically  derived  ground  state  potential  is 
reported  to  be  accurate  to  ±  0.05  eV  at  0.5  A  and  ±  0.01  eV  at  1.8  A.  The  third 
potential  is  that  derived  by  Leonas  and  Sermyagin  (Ref.  71)  from  differential 
elastic  scattering  cross  section  measurements.  For  internuclear  separations 
greater  than  1.8  A,  which  affect  the  normalization  of  the  continuum  nuclear 
wavefunctions,  we  used  an  accurate  ab  initio  potential  that  has  been  reported  by 
Liu  and  McLean  (Ref.  72).  This  latter  potential  was  smoothly  joined  to  the 
Foreman,  Ceperley  or  Leonas  potentials  in  the  region  of  1.6  -  1.8  A. 

We  have  evaluated  the  Einstein  A-coef ficients  for  He2  corresponding  to  the 
A1!^  ♦  X1!^  and  *  X1!*  transitions  for  the  wavelength  region  (60  nm  <  X  < 

110  nm) ,  which  covers  the  experimental  He2  continuum  emission  reported  by 
Huffman,  et  al.  (Ref.  60).  Using  the  A-coef ficients  as  a  direct  measure  of  the 
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observed  continuum  intensity,  we  show  in  Fig.  7  a  comparison  of  the  data  of 
Huffman,  et  al .  (Ref.  60)  with  our  calculated  intensity  profiles.  Our  calculated 
profiles  are  the  result  of  a  least-squares  fit  of  the  A1!*  +  X1!*  and  l)l£*  -*■  X*Eg 
A-coef f icients  assuming  contributions  from  the  first  four  vibrational  levels  of 
both  transitions. 

The  calculated  vibrational  population  of  the  A1!!*  and  Dl£*  states  are  given 
in  Table  II  for  the  three  different  ground  state  He2  potentials  that  were 
studied.  Overall,  both  the  Foreman  (Ref.  69)  and  Ceperley  (Ref.  70)  potentials 
exhibit  a  good  fit  to  the  data  of  Huffman,  et  al .  (Ref.  60).  The  potential  given 
by  Foreman,  et  al.  (Ref.  69)  yields  a  statistically  better  fit  with  y2  =  0.002  as 
compared  to  y2  “  0.01  for  the  Ceperley  and  Partridge  potential  (Ref.  70).  A 
significantly  poorer  fit  with  y2  =  0.47  is  found  for  the  potential  proposed  by 
Leonas,  et  al .  (Ref.  71).  All  other  ground  state  potentials  that  we  have 
explored,  including  those  reported  by  Jordan  and  Amdur  (Ref.  73),  yield 
significantly  poorer  statistical  fits.  We  have  not  analyzed  tne  new  potential  of 
Nitz,  et  al .  (Ref.  74),  which,  as  mentioned  above,  is  similar  to  that  of  Foreman, 
et  al .  (Ref.  69). 

The  best  fit  of  our  calculated  intensities  to  the  experimental  data  of  the 
Hopfield  helium  continuum  requires  a  repulsive  energy  of  3.4  ±  0.1  eV  for  the 
ground  state  of  He2  in  the  vicinity  of  1.04  A  internuclear  separation.  This 
result  is  based  on  an  average  of  the  Foreman  and  Ceperley  potentials  since  the 
experimental  uncertainties  in  the  measured  intensity  distributions  prevent  us 
from  making  a  statistically  favored  choice. 

On  the  basis  of  this  study  of  the  He2  continuum  spectra  reported  by  Huffman 
et  al .  (Ref.  60),  Michels  and  Peterson  concluded  that  the  Foreman  ground  state 
potential  for  He2  was  essentially  correct.  An  exhaustive  study  by  Michels  and 
Peterson  (Ref.  75)  of  possible  shapes  of  the  potential  energy  curves  for  He2“ 
resulted  in  the  conclusion  that  the  He2~  beam  must  have  been  rotationally  hot  in 
the  SRI  experiments.  A  suitable  fit  to  the  experimental  spectrum  was  finally 
obtained  by  assuming  a  rotational  temperature  of  ~  15000  K,  with  significant 
contributions  from  higher  vibrational  levels  of  He2~.  This  collaborative 
experimental  and  theoretical  effort  solved  a  formidable  problem  but  introduced 
new  questions  concerning  our  understanding  of  electron  autodetachment  processes. 
Our  conclusion  that  vibrational  autodetachment  must  be  slow  in  the  He2-  system 
runs  counter  to  previous  ideas  and  suggests  that  several  electron  coupling 
mechanisms  must  be  reexamined  theoretically. 


There  have  been  numerous  therraochemical  and  kinetics  studies  of  weakly  bound 
ionic  clusters  during  the  past  decade  (Refs.  76  and  77).  Of  particular  interest 


7VU 


R88-927258 

are  the  H*  cation  and  H“  anion  sequences  since  they  represent  examples  of  the 
chemistry  of  strong  chemical  bonding,  in  H^,  the  simplest  hydrogen  bonded 
molecular  ion,  in  H^,  and  the  simplest  atomic  anion,  H~.  In  addition,  for  the 
first  few  members  of  these  sequences,  detailed  ab  initio  quantum  mechanical 
calculations  of  their  equilibrium  bond  lengths  and  vibrational  spectra  can  be 
carried  out. 

Very  recently,  the  possibility  of  forming  large  H+  clusters  (n  >  1000),  that 
could  be  kept  for  long  periods  of  time  in  a  Penning  trap,  has  been  discussed 
(Ref.  78).  Such  a  system  is  a  prototype  for  bulk  antiproton  or  antihydrogen 
storage,  with  uses  as  high  energy  density  material  or  as  an  energy  source  for 
advanced  propulsion  concepts.  The  mechanism  for  formation  and  storage  of  these 
large  cluster  ions  can  be  understood  only  if  the  thermodynamics  and  kinetics  of 
the  cluster  growth  processes  are  identified  and  analyzed  in  detail.  The  purpose 
of  our  study  was  to  examine  the  thermodynamic  stability  of  the  first  few  members 
of  the  hydrogen  cation  and  anion  cluster  sequences  since  their  stability  is  a 
requirement  for  defining  a  useful  association  pathway  for  the  growth  of  large 
clusters . 

Numerous  experiments  (Refs.  79  and  80)  have  measured  the  stability  of  the 
odd-merabered  hydrogen  clusters  that  are  formed  by  the  association  reaction: 

Hn  +  H2  ^  Hn+2  ’>  n  =  3>5<7>--*  (H> 

These  clusters  exhibit  stabilization  energies  of  approximately  3-6  kcal/mol,  in 
good  agreement  with  theoretical  calculations  (Refs.  81  and  82).  Until  recently, 
the  even-membered  hydrogen  clusters  had  not  been  observed,  with  the  exception  of 
the  suggestion  of  a  stable  ion  formed  by  the  attachment  of  a  H  atom  to  the 

ion  (Ref.  83).  A  very  recent  study  of  the  even-membered  clusters  (Ref.  84) 
suggests  that  the  observation  could  be  interpreted  as  H2D+.  This  study  also 
indicates  a  high  degree  of  stability  for  . 

In  contrast,  the  negative  hydrogen  ion  clusters  have  never  been  observed  and 
a  recent  careful  search  for  stable  or  raetastable  and  produced  negative 
results  (Ref.  85).  There  have  been  several  theoretical  studies  (Refs.  86  and  87) 
which  support  the  existence  of  stable  HJ,  and  higher  clusters,  but  more  accurate 
ab  initio  calculations  (Ref.  88)  shed  doubt  on  this  conclusion. 

In  order  to  establish  the  thermodynamic  stability  of  these  hydrogen  cluster 
ions,  ab  initio  calculations  were  carried  out  to  determine  their  electronic 
structure  and  stable  geometries.  Several  Gaussian  basis  sets,  with  increasing 
flexibility,  were  employed  in  these  studies  and  levels  of  theory  up  to  single  and 
double  configuration  interaction  expansions  (CISD)  were  examined. 
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The  results  of  these  calculations  are  summarized  in  Table  III,  which 
includes  the  H"  and  H2  constituents  as  well  as  H^.  It  is  clear  from  Table  III 
that  the  inclusion  of  the  vibrational  zero  point  energy  is  necessary  to  assess 
the  stability  of  H^.  We  find  with  every  basis  and  level  of  theory  that  the 
ion  is  not  stable  with  respect  to  H-  and  H2 .  This  conclusion  is  supported  by  the 
results  of  earlier  theoretical  work  (Refs.  87  and  88),  and  by  the  recent 
experimental  study  of  Bae,  Coggiola,  and  Peterson  (Ref.  85). 


Our  calculations  appear  to  rule  out  the  formation  of  hydrogen  anion  clusters 
though  either  of  the  association  reactions: 


H~  +  H  +  H  7, 
n  n+ 1 


Hii  +  H2 


n+2 


(12) 

(13) 


Although  previous  studies  by  Hirao  and  Yamabe  (Ref.  88)  indicate  a  weak  stability 
for  H“  clusters,  up  to  n  =  13,  the  non-existence  of  and  precludes  formation 
of  higher  clusters  through  either  reaction  (12)  or  (13).  There  remains  the 
possibility  (albeit  remote)  that  association  of  H~  with  a  stable  Rydberg  state  of 
H3  could  produce  H^.  Further  calculations  are  in  progress  to  access  this 
association  pathway. 
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LOW-LYING 
Separated  atoms 
Li  +  Li 


TABLE  I 


MOLECULAR  STATES  OF  Li2  AND  THEIR  DISSOCIATION  LIMITS 


E(cm-1 ) 


—  1  * 


Molecular  states 


2V2s) 

+ 

2Sg(2s) 

0.0 

■r*g(D, 

3E*(1) 

2Sg(2s) 

+ 

2Pu(2p) 

14903.0 

lKU)> 

3yn. 

3E*(1) 

'V11- 

1nu(1), 

3VU- 

3nu(1) 

2Sg(2s) 

+ 

2Sg°S) 

27206.0 

1£*U). 

^u(i), 

3yn, 

3E^CD 

2Pu(2p) 

+ 

2pu(2p) 

29807.0 

^"(D, 

3z-(i) 

>IgU). 

lnud), 

3nga), 

3nu(i) 

Vl)- 

3au(D, 

2SgC2s> 

+ 

2pu(3p) 

30925.0 

1£g<1)- 

^u(1)« 

3e*U), 

3£u<1), 

V”* 

V1’- 

3JIU(1) 

2V2s) 

+ 

2Dg(3d) 

31283.0 

‘rjd). 

3lga), 

3y». 

llTg(1)« 

lnu!1)’ 

3ng(i), 

3nu(i). 

V1’* 

V11' 

VD, 

3Vl) 

Li2* 

+ 

e 

33200.0 

2V2s) 

+ 

2Sg(4s) 

35012.0 

l£sa)- 

^u(1)* 

3e*U). 

3e-(D 

2V2s) 

+ 

2Pu(4p) 

36469.0 

1 *;<!>. 

3yo. 

3r;a>, 

ing(i), 

■yn, 

3yi», 

X(1) 

C.  E.  Moore,  Atomic  Energy  Levels,  _1_,  9  (1949). 
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TABLE  1 

(Continued) 

(2s) 

+  ^  (4d) 

O 

36623.0 

^U), 

'V1'* 

\<‘l. 

3V°’ 

3nu(D, 

lv“’ 

V‘>- 

Vl>- 

2sg 

(2s) 

+  2Fu(4f) 

36630.0 

vu’ 

3Eg(1)' 

V1J. 

‘V11’ 

>nuto, 

3ngU). 
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■sua), 

3Vn' 

V1}- 

■Vn, 

3*gU), 

Vu 

2sg 

( 2s ) 

+  2S  (5s) 

© 

38299.0 

1£+g(1)> 

3E*o>, 

3e;u) 

Li+ 

Li~ 

[Js 
[  is 

(Is2)]  + 
g( 2s  2)  ] 

^474.0 

1Lg(1)> 

>i;o) 

2s8 

(2s) 

+  2Pu(5p) 

39015.0 

vn* 

V”. 

‘V‘>. 

ln„U>. 

3VU' 

3nu(1) 

2SS 

(2s) 

+  2Dg(5d) 

39094.0 

‘tJU). 

© 

Vl)» 

lVl>- 

Vu’ 

3”gll)' 

3nuu>. 

V‘>’ 

1au(D. 

3VU- 

VD 

'sg 

(2s) 

+  2Fu(5f) 

39104.0 

V1}’ 

3Z^(1), 

30‘>. 

Vl)* 

1nud). 

3ng(i), 

3Hu(1)' 

vn- 

Vl)* 

1V1)> 

‘♦u(U, 

3*g(l), 

3»„<d 

2sg 

(2s) 

+  Doublets 
(5<n?9) 

39987.0 

42003.0 

4 

Singlet , 

E ,  n,  L 

,  triplet 
.  (g,u) 
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2pu(2p)  +  2Sg(3s) 


2S  (2s)  +  Doublet 
(9<n<«) 


Li+  [IS  (Is2)]  + 
Li  [ 2Sg?2s) ] 


2Pu(2p)  +  2Pu(3p) 


2Pu(2p)  +  2Dg(3d) 


2Py(2p)  +  Doublet 


(n=4 ,  5) 


2S  (3s)  +  2S  (3s) 

O  O 


Li  +  PS  (Is2)]  + 
Li~  (  1Pu(2s2p) ] 


TABLE  1  (Concluded) 
42110.0 


42118.0 


43482.0 


45829.0 


46187.0 


51088.0 

53782.0 


54412.0 

57920.0 


30 


1Ej(l) ,  1Z^(1) ,  3E+(1), 


3z+(i),  ^gd),  1nu(i), 


3ng(i),  3nu(i) 


Singlet,  triplet 
E,  II,  4  (g,u) 


2Ej(l),  2Z^(1) 


xEg(2) ,  lE*(2),  3Eg(2) , 


3E£(2),  lE“(l),  1 E“( 1 ) , 


3E“(1),  3e"(D,  1ng(2) , 


1nu(2) ,  3ng(2) ,  3nu(2),  ^(l), 


1Au(l) ,  3Ag(l),  3Au(1) 


1Ej(2)  ,  1  E*(2) ,  iEgd), 


1 E~ (1 ) ,  3e:(2),  3r+(2), 


3Eg(l),  3e~(i),  1ng(3) , 


:HU(3),  3ng(3),  3nu(3) ,  1Ag(2), 


1AU(2),  3Ag(  2) ,  3  Au(2)  l*g(l), 


UyU).  3$g(l)>  3*U(l) 


Singlet,  triplet 
E,  n,  A,  *,  T  (g.u) 


li+d),  3E+d) 


^td),  iz+d),  Jn  (i), 

‘1*0) 
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TABLE  II 


Least-squares  Vibrational  Population  Analysis  of  the  Calculated 
Hopfield  Spectral  Distribution 


Vibrational  Level 


Ground  State  He„  Potential 


Foreman 

Ceperley 

Leonas 

< 

II 

o 

0.706 

0.739 

0.628 

V  =  1 

0.182 

0.166 

0.210 

V  =  2 

<0.002 

0. 

0.029 

v  =  3 

0. 

0. 

0. 133 

w; 


m. 


Table  III.  Calculated  Electronic  Energies  and  bond  Lengths  for 
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Spectral  intensity  distribution  of  the  Hopfield  helium  continuum. 
Comparison  of  experimental  and  calculated  data. 
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APPENDIX  A 

ELECTRONIC  STRUCTURE  OF  He”  * 

H.  H.  Michels 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 


ABSTRACT 

The  electronic  structure  of  the  He2  anion  has  been  studied  using  optimized 
Cl  wavefunctions .  We  find  that  the  ^11  state  of  He2  is  bound  relative  to  the 

a3!'*'  state  of  He.,,  with  a  calculated  electron  affinity  of  0.182  eV.  Only  the 

^  ...  .  r  3+i 

lowest  vibrational  level  of  this  anion  is  stable  relative  to  the  He2  [a  £uJ  +  e 

continuum. 


*Supported  in  part  by  AFOSR  under  Contract  F49620-85-C-0095 . 
Published  in  Chemical  Physics  Letters,  126,  1986,  pp.  537-540. 
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APPENDIX  B 


Spectral  Intensity  of  the  Hopfield  Helium  Continuum: 

Sir 

An  Analysis  of  the  Ground  State  Potential  for  He2  . 

H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 


J.  R.  Peterson 
SRI  International 
Menlo  Park,  CA  94025 


ABSTRACT 


The  intensity  distribution  of  the  Hopfield  helium  continuum  has  been 
analyzed  using  accurate  numerical  wavefunctions .  The  best  comparison  of  our 
calculated  intensity  to  the  experimental  continuum  of  Huffman,  et  al ,  yields  a 
repulsive  energy  of  3.4  ±0.1  eV  for  the  ground  X  state  of  He2  in  the 

vicinity  of  1.04  A  internuclear  separation. 


*Supported  in  part  by  AFOSR  under  Contracts  F49620-85-C-0095  and  F49620-85-K-001 7 
and  NSF  Grant  PHY-8410980. 

Published  in  Chemical  Physics  Letters,  134,  1987,  pp.  571-574. 
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APPENDIX  C 

«  "if 

The  Electronic  Structure  and  Stability  of  the  Anion 

H.  H.  Michels  and  J.  A.  Montgomery,  Jr. 

United  Technologies  Research  Center,  East  Hartford,  CT  06108,  USA 


ABSTRACT 

A  systematic  study  of  the  electronic  structure  of  the  anion  has  been 
carried  out  using  both  perturbation  theory  and  configuration  interaction. 
Counterpoise  corrections  are  used  to  estimate  the  basis  set  superposition  error 
in  the  computed  dissociation  energy.  We  find  that  is  not  thermodynamically 
stable  relative  to  H”  +  H2  upon  consideration  of  the  vibrational  zero-point 
contribution  to  the  dissociation  energy.  Isotopic  analysis  predicts  weak 
stability  for  and  D2-H“. 


♦Supported  in  part  by  AFOSR  under  Contract  F49620-85-C-0095 .  Published  in 
Chemical  Physics  Letters,  139 ,  1987,  pp.  535-539. 
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APPENDIX  D 

Negative  Ion  Formation  in  Lithium  Atom  Collisions 

H.  H.  Michels  and 
J.  A.  Montgomery,  Jr. 

United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 

ABSTRACT 

The  formation  of  Li-  by  dissociative  attachment  in  e  +  Li2  collisions  is 
characterized  by  a  large  cross  section  for  electron  attachment  to  highly  vibra- 
tionally  excited  Li2  molecules.  However,  the  electronic  structure  of  the  Li2“ 
system  dictates  that  low  energy  electrons  (<  0.2  eV)  are  needed  to  minimize  loss 
processes.  An  alternate  production  route  for  Li-  is  possible  via  direct  ion-pair 
formation  through  collisions  of  highly  excited  (Rydberg)  Li  atoms.  The 
mechanisms  for  ion-pair  formation  will  be  discussed  in  terms  of  the  high-lying 
molecular  Rydberg  states  of  Li^. 


★Supported  in  part  by  AFOSR  under  Contract  F49620-85-C-0095 . 

Published  in  the  Proceedings  of  the  Fourth  International  Symposium  on  the 
Production  and  Neutralization  of  Negative  Ions  and  Beams  (Brookhaven,  1986),  ed. 
J.  Alessi,  1987. 
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APPENDIX  E 

Electronic  Structure  and  Stability  of  Small  Cation 
and  Anion  Hydrogen  Clusters 


■jlf 

J.  A.  Montgomery,  Jr.  and  H.  H.  Michels 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 

ABSTRACT 


Ab  initio  calculations  of  the  electronic  structure  of  H  ”  and  H  +  clus- 
-  n  n 

Lets  nave  been  carried  out  using  accurate  Gaussian  basis  sets  and  with  levels 
of  theory  up  to  fourth-order  perturbation  theory  (MP4)  and  single  and  double 
excitation  configuration  interaction  (CISD).  The  odd  hydrogen  cation  addition 
sequence,  Hn+  +  H2  -*■  Hn*2,  appears  to  be  thermodynamically  stable  for  large 
size  cluster  formation.  The  even  Hg+  cation  also  exhibits  surprising  stabil¬ 
ity  in  D2c(  symmetry.  In  contrast,  the  hydrogen  anion  addition  sequences,  Hn" 

+  H  -*■  Hn“j  and  Hn"  +  H2  -»  Hn“2 ,  aPPear  to  be  thermoneutral  or  unstable. 


♦Supported  in  part  by  AFOSR  under  Contract  F49620-85-C-0095  and  in  part  by 
AFRPL  under  Contract  F04611-86-C-0071 . 


Published  in  the  Proceedings  of  the  1987  Cooling,  Condensation  and  Storage 
of  Hydrogen  Cluster  Ions  Workshop,  J.  T.  Bahns,  ed . ,  1988,  pp.  83-94. 
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APPENDIX  F 

**e2  [** Ilg ]  +  He2  [X^Eg+]  Autodetachment 

Energy  Spectrum:  Assessment  of  the  He2~  and 
He2  Ground  State  Potentials 

Y.  K.  Bae  and  J.  R.  Peterson 
SRI  International 
Menlo  Park,  California  94025 


H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  Connecticut  06108 


ABSTRACT 


The  energy  spectrum  from  the  spin-forbidden  electronic  autodetachraent  from 
the  metastable  He2  [^11  ]  ion  has  been  measured.  The  single-peaked  continuum, 
resulting  from  a  vertical  transition  to  the  repulsive  wall  of  the  He2  [x1!*] 
ground  state,  has  a  maximum  at  15.70  ±  0.15  eV.  For  a  rotationally  cool  beam, 
the  spectrum  is  inconsistent  with  the  existing  potentials  of  He2~  and  He2  [X1!*]. 
A  suitable  fit,  neglecting  rotational  effects,  would  require  an  increase  of 
0.06  A  in  the  calculated  Rg  of  He2  •  However,  a  very  good  fit  can  also  be 
obtained  from  the  existing  potentials  by  assuming  a  rotationally  hot  (15000  K) 
distribution.  Surprisingly,  both  fits  require  contributions  from  levels  v  >  0, 
which  are  known  to  undergo  vibrational  autodetachment  to  He2[a"^£*].  For  such 
contributions,  the  latter  process  must  be  exceptionally  slow  (<10^/s). 


Supported  in  part  by  NSF  under  Grant  PHY-84-10980  and  AFOSR  under  Contracts 
F496 20-85 -C-00 17  and  F49620-85-C-0095 . 

Published  in  Physical  Review  A,  37^  1988,  pp.  2778-2784. 
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